Electron transport in a quantum wire with realistic Coulomb interaction 
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Electron transport in a quantum wire with leads is investigated with actual Coulomb interaction 
taken into account. The latter includes both the direct interaction of electrons with each other 
and their interaction via the image charges induced in the leads. Exact analytical solution of the 
problem is found with the use of the bosonization technique for one-dimensional electrons and the 
three-dimensional Poisson equation for the electric field. The Coulomb interaction is shown to change 
significantly the electron density distribution along the wire as compared with the Luttinger-liquid 
model with short-range interactions. In dc and low-frequency regimes, the Coulomb interaction 
causes the charge density to increase strongly in the vicinity of the contacts with the leads. The 
quantum wire impedance shows an oscillating behavior versus the frequency caused by the resonances 
of the charge waves. The Coulomb interaction produces a frequency- dependent renormalization of 
the charge wave velocity. 

PACS number(s): 72.10.Bg 
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I. INTRODUCTION 



The electron-electron interaction is generally rec- 
ognized now to be fundamentally important in one- 
dimensional (ID) structures. The interaction of ID elec- 
trons turns out to be so significant that the Fermi liquid 
concept breaks down. More adequate becomes the no- 
tion of a strongly correlated state known as Luttinger 
liquid (LL) with bosonlike excitationsIM The commonly 
used LL model treats the electron-electron interaction 
as the short-range one. In this approach, the electron- 
electron interaction changes the electron liquid compress- 
ibility giving rise to a renormalization of charge- wave ve- 
locity. Generally speaking, the renormalization param- 
eter g is a function of the wave vector p of boson ex- 
citations. However, within the short-range interaction 
approach the parameter g is supposed to be a constant. 
Presently there is no unambiguous evidence that the LL 
really exists in quantum wires or those rejecting this con- 
cept. The experiment of Tarucha et alB has shown that 
the dc conductance of a quantum wire structure is quan- 
tized by standard steps e 2 /h. JJtis result was explained 
in thc.£rame of the LL model.aU More recently, Yacoby 
et al.ua have found a nonuniversal conductance quanti- 
zation by steps different from e 2 /h. One can therefore 
conclude that the experiments on quantum wires reveal 
more complex behavior of ID conduction than the sim- 
ple LL model predicts. Thus an important problem is to 
develop the theory for actual quantum wire structures. 
Application of the LL theory for this purpose points out 
two problems. 

First, the assumption that electrons interact with each 
other locally is evidently inadequate in the real situa- 
tion since the Coulomb interaction is-essentially nonlo- 
cal. This assumption is often justifiectalllJ by the screen- 



ing effect of the highly conducting gate electrode. la 
this case the gate current should be taken into account!!!! 
to provide the charge ccinservation and cause the the- 
ory to be gauge invariant £3 Hence the screening effect of 
conducting electrodes (gates and leads) should be thor- 
oughly analyzed and taken into consideration in order 
to understand the experimental situation. In reality, the 
screening by the electrodes consists in the appearance 
of image charges of the electrons which are situated in- 
side the quantum wire. Due to the image charges, the 
electron-electron interaction becomes dipolelike (or gen- 
erally multipolclikc) but the dependence of g on p may 
be essential for the structures realized experimentally. 

Second, the conductanij£-.of the mesoscopic quantum 
wire structure is knownotJ to be substantially deter- 
mined by the contacts of the wire with the leads. Besides, 
the contacts between the quantum wire and the leads 
cause reflection of bosonlike excitations which deter-minps. 
the high-frequency behavior of the admittance.Ba'E£l 
Thus, the interaction of electrons moving in the quan- 
tum wire with leads has to be taken into consideration. 

The purpose of the paper is to obtain the actual form 
of the electron-electron interaction potential in quantum 
wire structures with leads and to investigate the phase- 
coherent transport of electrons in both the dc and ac 
regimes. The difficulty of this problem is caused by a 
nonlocal nature of the interaction and by the fact that in 
a quantum wire of finite length the translational symme- 
try is broken and hence the electron-electron interaction 
potential depends separately on the coordinates x and 
x' of interacting electrons rather than on the difference 
|x — x'\. We have found a situation in which this prob- 
lem is solved exactly in the frame of the bosonization 
technique. It is realized when the lead surfaces may be 
approximated by planes perpendicular to the wire. 

It is worth noting that there is an alternative (con- 
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tactless) approach in studying the ac transport of elec- 
trons in quantum wires. In this case the quantum wire 
is not supplied with leads. The ac transport is investi- 
gated by means of measuring the absorption or scattering 
of the electromagnetic radiation. This situation was r£=. 
cently considered by Cuniberti, Sassetti, and KramercEl 
for a homogeneous quantum wire of infinite length. They 
have investigated ac conductance defined via the absorp- 
tion of electromagnetic radiation taking into account the 
electron-electron interaction of finite range and an arbi- 
trary distribution of the external electric field along the 
wire. It was found that both the interaction length and 
the electric field distribution affect significantly the ac 
conductance. In the present paper we show that the leads 
produce an essential effect due to inhomogeneity of the 
electron-electron interaction in the wire. It manifests it- 
self in the charge-density distribution along the wire and 
the frequency dependence of the impedance. 

The paper is organized as follows. In Sec. II the po- 
tential of electron-electron interaction in a quantum wire 
structure with leads is obtained for the electron inter- 
action with the charges induced on the leads taken into 
account. Section III describes the equation of motion 
for bosonized phase field with nonlocal interaction and 
gives its solution via expansion in terms of the eigen- 
functions of the electron-electron interaction potential. 
In Sec. IV, charge-density distribution in the structure 
is investigated for both the dc and ac regimes. Section 
V contains the calculation and analysis of the impedance 
of the quantum wire structure. 

II. ELECTRON-ELECTRON INTERACTION 
POTENTIAL 

The mesoscopic structure under consideration consists 
of a quantum wire coupled to two bulky (2D or 3D) re- 
gions (the electron reservoirs) which serve as leads. The 
electrons in the wire interact with each other both di- 
rectly and via the surface charges which are induced on 
the surface of the leads. The electron-electron interaction 
energy W is defined by the product of the electron den- 
sity p(r) at a point r and the potential i/3,(r, r') created 
at this point by the charge at a point r'. This poten- 
tial is determined by the Laplace equation with bound- 
ary conditions corresponding to the given configuration 
of the leads. When calculating <^(r,r'), it is reasonable 
to consider the lead surfaces as equipotential ones. This 
is a natural assumption. As we are interested mainly in 
the electron behavior in the quantum wire, we can as- 
sume that the characteristic times of electron processes 



(such as Maxwell relaxation and plasma waves) inside 
the reservoirs are much shorter than the electron transit 
time through the quantum wire. This will be the case if 
the reservoirs are perfectly conducting. 

Distribution of the electron density p(r) in the channel 
can be written in the form 

p(r) = x(r±)p(x) , (1) 

where x( r J-) is a normalized function of the radial coordi- 
nate perpendicular to the channel and p{x) is a function 
of the coordinate along the channel. 

The potential ipi(r, r') created by an electron density is 
expressed via the Green function of the Laplace equation 
with zero boundary conditions at the lead surfaces. 

The product of p(r) and fi(r, r') can be integrated 
over the transverse coordinates and r'j_ to give the 
following expression for the electron-electron interaction 
energy O 

W = i J I dxdx' p(x)p(x')U(x,x'). (2) 
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FIG. 1. Schematic view of a quantum wire structure with 
leads, a moving charge, and electric force lines. 

An explicit form of U(x,x') was found in Ref. ^ for 
realistic situation where the electrode surfaces are two 
planes x = —L/2 and x = L/2 perpendicular to the 
channel (Fig. 1). In this case 

U(x,x') = - [ d 2 q\ Xq \ 2 G q (x,x'), (3) 
ne J 

e is the dielectric constant outside of the channel, Xq is 
the Fourier-Bessel transform of the radial function x( r ±)i 
and 



Gq{x, X ) 



L 



qLsmh(qL) 



smh[q(L/2 + x)] sinh[q(L/2- x')} if x < x' 
smh[q(L/2-x)}sinh[q(L/2 + x')] if x > x' 



(4) 
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The interaction potential denned by Eqs. (|^) and (0) is 
shown in Fig. 2 as a function of x for a variety of x' , with 
V(x,x') being a normalized form of U(x,x'): 

V(x,x') = U(x,x')eL/e 2 . 



cific configuration of the leads. The potential defined by 
Eqs. (||) and (||) will be used below in getting an exactly 
solvable model of the interacting electron transport. 



III. THE EQUATION OF MOTION 
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FIG. 2. Distance dependence of the electron-electron inter- 
action potential for a variety of x'. Dotted line is the envelope 
of V(x, x) maxima. 

When the intcrclectron distance is larger than the 
wire radius a (|ar — x'\ ^> a), the potential decreases as 
V ~ L/\x — x'\. In the middle part of the quantum wire 
V(x = x') ~ L/a. Near the contacts (x,x' — > ±L/2), V 
goes to zero due to screening effect of the charges induced 
on the lead surfaces. The behavior of this kind is quite 
general for the interaction potential regardless of the spe- 



To study a linear response of the quantum wire struc- 
ture to an external voltage V a exp(-iujt) applied across 
the leads, we can restrict ourselves to the consideration 
of low-energy excitations of the electron system. The 
most adequate methodriar this purpose is the standard 
bosonization technique .Era We will use this technique as- 
suming that the electron density fluctuations are long- 
range ones. 

When external voltage is applied across the electrodes, 
the electric potential p(r) in the wire is determined by 
the Poisson equation with the boundary conditions con- 
trolled by the applied voltage: <p — at the left reservoir 
and tp — V a at the right reservoir. It is convenient to 
present the electric potential in the wire as a sum 

<P = Pcxt + (pi ■ (5) 

Here ip ex t is the potential which would be without the 
wire. It is determined by the Laplace equation with the 
same boundary conditions as the total potential ip. The 
potential ipi is defined by the Poisson equation with zero 
boundary conditions. The potential ipi is precisely the 
one used in Sec. II while calculating the electron-electron 
interaction energy [see Eqs. (||) and 

Thus the bosonized HamiltonianffEa of ID spinless 
electrons can be taken in the form 



H = 



' dx 
2^ 



(l+.gi) 7 r 2 n 2 + (l- 3l )(9 :E $) 2 -e / dx p(x, t)ip ext (x, t) + 



dx dx' 
2tt 2 



{d x $)U(x,x')(d x ,$), (6) 



where t) is the phase field related to the charge excitations, II(a;, t) is the momentum density conjugate to vf 
is the Fermi velocity, and g\ is the backscattering parameter.tSI _ 

By writing this Hamiltonian, we assume implicitly that the ground state is uniform. More careful investigationlij 
shows that really the ground state is nonuniform due to two factors: (i) charging of the quantum wire which occurs 
as a consequence of the electron transfer between the wire and the reservoirs during the process of the establishment 
of the equilibrium electrochemical potential and (ii) Friedel oscillations near the contacts. The charge stored in the 
wire depends on the wire radius and the background density of the positive charge. This charge may be both positive 
and negative. Under certain conditions the wire remains neutral. The charging effect will be investigated elsewhere. 
In the present paper the ground state is assumed to be neutral. The Friedel oscillations of the electron density have 
a characteristic length of the order of the Fermi wavelength. Since we restrict our consideration by the long-range 
variation of the electron density, the ground state may be considered as uniform. 

The long wave component of the electron density p is related to $ by 

p=~d x $. 

From the Hamiltonian (ph we get the following equation of motion for the phase field $: 



L/2 



e p C xt + - 

7T 



dx'U(x,x')d x ><$> 



-L/2 



(7) 
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with 



fj = 



.91 



l-5i 



Following Refs. 4-6, we extend the one-dimensional 
Eq. (0) to the reservoirs assuming that the electron den- 
sity inside them is extremely high and their conductivity 
is ideal. In such a case, tp cxt is a constant inside the 
reservoirs. Hence the first term on the right-hand side 
of Eq. (0) can be omitted. The second term (appear- 
ing from the electron-electron interaction energy W) is 
known to be small as compared with the kinetic energy 
on the left-hand side of Eq. (^) when the electron density 
is high. In this sense, the electrons in the reservoirs are 
nonintcracting although the external field is of course ide- 
ally screened there. Thus, the right-hand side of Eq. ([?]) 
may be dropped in the reservoirs. The solution of the 
bosonized equation in the reservoirs should satisfy the 
condition that the density wave be restricted at |x| — > oo 
when the external voltage is turned on adiabatically. The 
boundary conditions at the contacts between the wire 
and the leads require the continuity of $ and the par- 
ticle current. The latter means that (v/g)d x $ must be 
continuous there. 



We find the exact solution of Eq. (Q) with the electron- 
electron interaction potential U (x, x') defined by Eqs. (||) 
and (0). In terms of dimensionless variables 



v $ 
egLV a 



f = 



v a 



Eq. (Q) takes the form: 



d 



du » du 



n 2 u = o . 



(8) 



Here 



= 



1 ~ 

eg 



n = 



luL 



the operator V is defined as 

1/2 
-1/2 

where 



(9) 



dy 
smhy 



\X V \ 2 



sinh[y(l/2 + 0]sinh[y(l/2-r)] if £ < C, 
sinh[y(l/2-£)] sinh[y(l/2 + £')] if £ > 



(10) 



by 



The electron density in standard units is related to 



p(0 



2egV a 



H(0 



(11) 



After solving Eq. (g) in the reservoirs we use the conti- 
nuity of the electron flow and the phase <E> at the contacts 
to obtain the boundary conditions directly to Eq. (||) in 
the inner region of the quantum wire, —1/2 < £ < 1/2, 



± i fig itL 



Tl/2 







(12) 



The integro-differential equation of the form Eq. (|8j) 
may be solved via an expansion in terms of the eigen- 
functions of the operator V. It is easy to verify that the 
functions 



ipn(0 = v^sin 



are the eigenfunctions of V defined by Eqs 
with the eigenvalues 



1,2,3,. 
(0) and (II 



Ar. 



ydy 



\Xy\ 



(13) 



la V 2 + ( n7r ) 2 

For the geometry of the sample under consideration, 
the external potential /(£) is a linear function which may 
be expanded in terms of ip n (Q with even n. 



The exact solution of Eq. (|8|) with the boundary con- 
ditions ( |l2| ) can be obtained via the expansion in terms 
of ip n (£)- As a result, we get the following expressions 
for the dimensionless electron density ««;(£) : 

oo 

u dO = J2 c « sin + V 2 )] (14) 

n=l 

and the phase field u(£), 

oo 

u(0 = A(Sl) - ^ cos [2*n(£ + 1/2)] . (15) 

n— 1 

where 



c n = B(Q) 



Ann 



47r 2 n 2 (l + /3A 2 „)-0 2 



m= i i-«gHD(n) 



D(Sl) = J2 



n 2ig + n- 4ign 2 D(n) 

-2ig 

2ig + fl- AigVL 2 D{VL) ' 
1 



-47r 2 n 2 (l + /3A 2 „)-f2 2 



(16) 
(17) 
(18) 

(19) 
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IV. ELECTRON DENSITY DISTRIBUTION 



contact as smh(x/£), with characteristic length being 



The purpose of this section is to clarify how the real 
Coulomb interaction affects the value and the distribu- 
tion of the charge density in the quantum wire structure. 

First of all, let us consider a limiting case of the short- 
range interaction. In this case, V(x, x') oc S(x — x') and 
hence the eigenvalues are independent of n and A2 n = A. 
All the sums are easily calculated, which results in the 
following expression for the normalized electron density: 



«c(0 



1 



g* sin(fi*0 



Qy/TTM g* cos(r2*/2) - i sin(fi*/ 2 ) ' 



(20) 



where renormalized values g* — g/y/1 + (3X and f2* — 
Q/s/l + (5\ are introduced. The density /?(£) calculated 
according to Eq. (|2(]) coincides exactly with that found 
in Ref. [if] in the framework of the standard LL model 
with the interaction parameter g = g* ■ 
In the limit of Q -> 0, Eq. (poh yields 



h(0 



1 + /3A 



(21) 



With increasing, frequency the charge waves appear 
which have resonancesE3 along the wire when = 
2nn^/T+JX. 

Another case will be useful in what follows as a ref- 
erence point to demonstrate the Coulomb interaction ef- 
fect. It is the case of noninteracting electrons which cor- 
responds to (3 = in Eqs. @ and @. 

For the case of the realistic interaction, the electron 
density distribution given by Eq. ( |l4| ) is generally more 
complicated. However, simple results are obtained for 
the regions near the contacts and in the middle part oi 
the wire, taking into account the specific behavior of A„ 
versus n. It is determined by the fact that the radial 
function x( r -0 is located in the region of radius a which 
is much shorter than the wire length L, i.e., a = a/L 
is a small parameter. The results which will be given 
below are qualitatively valid for any localized function 
x(r_i_). To be specific, we will use the Gaussian form for 
x( r x) when it is necessary to bring the calculations to 
final form. One obtains that A„ varies slowly with n for 
nan <§C 1 and A„ decreases as n~ 2 for nan ^> 1. 

In the vicinity of the contacts, the main contribution 
to the sum in Eq. ( |l4| ) is due to the large-n terms for 
which an asymptotic expression of A„ ~ n~ 2 can be used. 
Thus the following expression for the normalized electron 
density is obtained in the vicinity of the left electrode 
[(1/2 + «1]: 



H(0 



sinh(£v/2/3/a 2 - ft 2 ) 
~ sinh(y2/3/a 2 -^ 2 /2) 



(22) 



with S(O) being defined by Eq. (|l|). Equation @ 
shows that (i) U£ decreases with the distance from the 



y / 2/3~a 2 n 2 ' 

(ii) at 17 = 0, the boundary value ofyUg is equal to 1/2 
and is independent of the interactionO 

In the middle part of the wire, Eq. ( |l4| ) may be sim- 
plified when a is exponentially small, i.e., -lna > 1. In 
this case the sum in Eq. ( |l4| ) may by estimated assuming 

A„« A = -ln(2^ 2 a 2 )- 7 , 

7 « 0.57721 . . . being Euler's constant. This calcula- 
tion results in the same equation as Eq. ( |20| ) for the 
short-range interacting electron gas, where A should be 
replaced by Ao- This is equivalent to introducing an ef- 
fective interaction parameter 



ffoff 



Vi + Wo 



into the LL model. 

For the sake of simplicity, we suppose g 



1 below. 
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FIG. 3. Distribution of the dc electron density p (thick 
solid line) and the potential ip (thin solid line) along the wire 
for long-range Coulomb interaction. The dashed line is the 
density distribution according to the LL model with g e ff, the 
dotted line is that for noninteracting electrons. The param- 
eters used are: the Fermi energy ef = 5 meV and a = 0.02. 
The interaction parameter /3 was calculated to be 0.35. 

The effect of the Coulomb interaction on the electron 
density distribution in the quantum wire is demonstrated 
by Fig. 3 for the dc condition when a positive potential 
is applied to the left electrode with respect to the right 
one. Under the action of the external electric field, the 
electron system in the quantum wire is polarized: the 
electron liquid is compressed in the left part of the wire 
and decompressed in the right one. Respectively, at the 
left end of the quantum wire the excessive electrons ap- 
pear while at the right end the electron density is de- 
creased. If the electron-electron interaction is omitted 
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[see the dotted line in Fig. 3 and Eq. (|l|) with = 0], 
the electron density decreases linearly with the distance, 
with the boundary value of the normalized density being 
equal to ±1/2. When the short-range interaction with 
the effective interaction parameter g e g is turned on, the 
density distribution remains linear (dashed line) . But the 
slope goes down as a result of a decrease in the electron 
liquid compressibility. Note that the boundary value ol 
the electron density also decreases. 

When the actual Coulomb interaction is turned on, 
the electron density distribution is changed qualitatively 
(thick line in Fig. 3). As compared with the noninter- 
acting case, the charge density in the middle part of the 
wire decreases, which may be interpreted as neutraliza- 
tion of the negative and positive charges due to their 
mutual attraction. Towards the contacts the charge den- 
sity increases reaching ±1/2 at the boundaries. This be- 
havior may be understood from the fact that near the 
contacts the charges in the wire are neutralized by the 
image charges in the electrodes. 

On the other hand, if we compare the Coulomb inter- 
action case with the short-range interaction model, we 
find that the actual Coulomb interaction leads to an in- 
crease of the electron density fluctuation near the con- 
tacts. This fact may be interpreted as a result of the 
decrease in the interaction parameter g due to screening 
the electron-electron interaction by the electrode. 

At finite frequency, this near-contact effect of the 
Coulomb interaction is preserved up to a characteristic 
frequency Q w = \/2/3 /a, above which the exponentially 
decreasing part of p(x) disappears. 

The electric potential ip is easily calculated using 
Eq. (||), Green function (|]), and the electron density p(r) 
found above. Distribution of ip(£) along the quantum 
wire is shown in Fig. 3. A good proportionality is found 
between p and ip for dc conditions when a/ L <C 1: 



« const -I eip(^) . 

7TV 



The main effect observed when the frequency grows is 
the appearance of the traveling charge waves. This is il- 
lustrated by Fig. 4, where the real part of the normalized 
electron density is shown as a function of the distance 
for a number of frequencies. The curves shown in—Fig. 4 
were obtained by numerical calculation of Eq. (pf) The 
frequency is given in a normalized form 

luL 

V = 

llTV F 

as labels to the curves. 
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FIG. 4. Electron density distribution along the wire for a 
set of frequencies. The curves are labeled by the normalized 
frequency v. The parameters used in the calculations are the 
same as in Fig. 3. 

Our analytical solution given by Eqs. Q and @ 
shows that there is a set of characteristic frequencies Q n 
which are determined by the poles of £>(£!): 

Q n = 27rnV 1 + PM n ■ (23) 

For fl = £l n , the electron density is a standing wave with 
zeros at the contacts, n being the wave number: 

fdu\ sin[27rn(£+ 1/2)] 
\d{) n K 4nn{l + /3\ 2n ) ' 

It is noteworthy that the electron flow turns to zero at the 
contacts simultaneously with the electron density. Thus 
O n is the resonant frequencies of the charge waves in the 
quantum wire. Under the resonant condition the electron 
density perturbation is locked inside the wire. 

Equation ( ^3| ) shows that the frequency O n depends 
on the wave number 2im/L in a nonlinear manner due 
to the dependence of A ra upon n. For n — ► oo, the reso- 
nant frequency is proportional to the wave number which 
corresponds to the soundlike dispersion. For low wave 
numbers [n < (na)^ 1 ], fl n is noticeably higher than one 
expects from the soundlike dispersion. 

It is instructive to compare Eq. ( p3| ) with the disper- 
sion equation for charge waves-4n an infinite quantum 
wire which was found by Schulzta using the bosonization 
technique, 

fi(p) = P VI + P V p > ( 24 ) 

where p is the wave vector normalized by L~ x and V p is 
the Fourier transform of the interaction potential, which 
is approximated by the modified Bessel function Ko(ap). 
A similar dispersion law was obtained by Das Sarma and 
HwangEil for ID plasmons in the long-wavelength limit 
in the frame of the random-phase approximation taking 



G 



into consideration the more general form of the interac- 
tion potential. Equation (|24|) differs from Eq. (|23| ) in the 
replacement \2 n by V p . A reasonable approximation for 
Mn is 



X 2n w exp (27r 2 a 2 n 2 ) £a (2?r 



2 2 2\ 

arn 



(25) 



where £i (z) is the exponential integral. It is worth not- 
ing that in the limit L — > oo, Eq. (25) is the same as 
the Fourier transform of the interaction potential used in 
Rcf. |l^ when the screening length is much larger than the 
quantum wire diameter. On setting p = 2irn and com- 
paring the expressions (|23| ) and (|3|), one can see that 
they are close when iran <C 1 and differ significantly in 
the opposite case. 

One can say that Eq. ( p3| ) is a discrete version of the 
dispersion relation for ID electrons which takes correctly 
into account the electron-electron interaction in a finite 
ID system. It will be shown below that the discrete char- 
acter of the resonant frequencies f2 n of finite quantum 
wire results in strong peculiarities of the frequency de- 
pendence of admittance. 



In what follows, the impedance is analyzed rather than 
the admittance which is usually considered, because the 
frequency dependence of the impedance shows more pro- 
nounced features caused by the charge waves. The real 
part of Z is 



p r7 _ h 1 

e 2 1 + [AgnD(Q)] 2 



(28) 



When O — ^ 0, the impedance is equal to h/e 2 and is 
independent of the interaction parameter 0. The fre- 
quency dependence of ReZ is mainly governed by that of 
D(Q). The resonant frequencies tt n are by definition the 
poles of D(n). Between the neighboring poles, there is a 
zero of D(n). Equation (p|) shows that ReZ = when 
\D\ -> oo and ReZ = h/e 2 when D -» 0. Thus with 
increasing frequency, ReZ oscillates between zero (which 
occurs at the resonant frequencies) and h/e 2 , these limit- 
ing values being independent of the interaction. The fre- 
quency dependence of ReZ is illustrated by Fig. 5, where 
three cases are compared: noninteracting electrons, the 
LL model with short-range interaction, and the electrons 
with actual Coulomb interaction. 



V. THE IMPEDANCE 



The electron current in a quantum wke is determined 
by the time derivative of the phase field.El In the terms oi 
the normalized phase u, the current is 



., . 2iu;e 2 gL . 

]{x,u)) = u(x,uj)V a . 

hv 



(26) 



The current calculated in such a way depends on the 
coordinate x along the wire. However, the electric cur- 
rent jmeas which is detected by a measuring device is 
obviously independent of x. This current is defined as a 
charge flow through the leads. Its formation process is 
analyzed in Appendix A as applied to the specific quan- 
tum wire structure considered here. The measured cur- 
rent is the sum of a trivial capacitance current and the 
current caused by the quantum wire presence. According 
to the Shockley theorem,! 2 ^! the latter is 



1 



L/2 



dx j(x, u>) . 



-L/2 



Due to Eqs. flU) and @), the current j u becomes 
e 2 

j u =i—2gA(n)V a , 

with A(fl) being defined by Eq. (Jl7|). This results in 
the following expression for the quantum wire structure 
impedance: 



in 



l-4ignD(n) 2g 



(27) 




2 3 

norm, frequency, v 

FIG. 5. Real part of the impedance versus the normalized 
frequency v for true long-range interaction ( — ), LL approach 
with g=0.63 ( ), and noninteracting electrons (••■)• 

An interesting result is that ReZ = for the reso- 
nant frequencies. Under resonant condition the time- 
dependent evolution of the electron density is essentially 
oscillations between the two ends of the wire, and elec- 
trons are not emitted/absorbed by the contacts. As 
a consequence, the electric current component in-phase 
with the applied voltage vanishes, and the real part of 
both the admittance and the impedance turns to zero. 
However, in reality there is a finite dissipation which was 
not taken into account. The inclusion of the dissipation 
into calculations should restrict the minimum of ReZ by 
some low value. 

The resonances of ReZ occur when the frequency is 
a multiple of the inverse time of flight of electron ex- 
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citations along the quantum wire. This-conclusion was 
confirmed for both the noninteractingCTLJ electrons and, 
the short-range interacting electrons in the LL model.ta 
The fact that in the case of the Coulomb interaction 
the impedance oscillations are nonperiodic may be inter- 
preted as a result of the frequency-dependent renormal- 
ization of the charge-wave velocity due to the Coulomb 
interaction. At low frequency, the velocity renormalized 
by the Coulomb interaction is essentially larger than vf- 
The resonance frequency spectrum shows that the veloc- 
ity decreases with frequency. It is worthwhile to note 
that the phase velocity is important since the resonant 
conditions are obviously related to the wave interference. 

The imaginary part of the impedance may be presented 
in the form 



ImZ = —ujL c 



h AgVLD(VL) 

e* 1 + [4ijfiD(ft)] 2 ' 



(29) 



where 



hL 

2e 2 v F 



N 




2 3 4 5 

norm, frequency, v 

FIG. 6. Imaginary part of the impedance versus the nor- 
malized frequency v. 

The frequency dependence of ImZ consists in the linear 
decrease caused by the first term in the brackets and os- 
cillations around this dependence due to the second term. 
This behavior is illustrated by Fig. 6. The linear depen- 
dence of the ImZ on the frequency is obviously dominat- 
ing, which allows one to-|interpret L c g as a frequency- 
independent inductance.113 

When the frequency is small, the second term on the 
right-hand side of Eq. ( p9[ ) is comparable with the first 
one. In this case we can expand Eq. (E9h , 



ImZ 



h CI 
^Tg 



-1 



2.9 



_ 2 oo 



E 



i 



n 2 ^ n 2 (1 + /3A2n) 



to estimate the second term. If the electron-electron in- 
teraction is omitted (f3 = 0, g = 1), the second term is 



equal to 1/3. Upon increasing j3, the second term de- 
creases and ImZ remains negative in spite of the fact 
that the backscattering parameter g slightly increases 
when the interaction is turned on. Thus the reactive 
part of the impedance is always inductive if the electron- 
electron interaction is repulsive independent of the inter- 
action strength. 

The behavior of the impedance we have obtained here 
correlates with that found for the case of both the 
short-range interactionllj and noninteracting electrons .Eil 
Somewhat different behavior of the impedance was found 
recently in Ref. [l^ for a quantum wire without leads 
based on a rather general approach, which allows one to 
consider an arbitrary distribution of the external electric 
field along the quantum wire. In this case the impedance 
was shown to include both the inductive and capacitive 
components. This difference originates from the different 
experimental situation which was considered. In Ref. |l6| , 
a homogeneous wire of infinite length with a continuous 
spectrum of eigensolutions was examined which results 
in a dispersion relation u = u>(k). The resonant fea- 
ture of the impedance is caused by the inflection point 
of u>(k) where the group velocity reaches a minimum, 
with the characteristic wavelength of the charge waves 
being of the order of the wire radius. In the present pa- 
per, we consider a more specific situation of a finite wire 
restricted by leads with a discrete spectrum of eigenfunc- 
tions. The resonances we have found are attributed to 
the finite length of the wire. They appear when the char- 
acteristic wavelength of the charge waves is of the order 
of the wire length, i.e., the frequency is much lower than 



the resonant frequency which appears in Ref. 16 



The total admittance of the quantum wire structure is 
formed by both the quantum wire impedance Z defined 
by Eq. (E7|) and the interelectrode capa citan ce Co which 
is necessarily present there. Using Eq. (A.S) one obtains 



v A 1 ■ r RcZ-i(lmZ + LoC \Zf) 
Y tot (Lu) = — - iu)C = — 12 . (30) 



It is of interest to find the eigenfrequencies of the ad- 
mittance (the impedance) which describe the behavior 
of the system under consideration as an element of elec- 
tric circuit. They are known to be determined by the 
poles and zeros of the admittance (the impedance). The 
admittance zeros characterize the system when the ex- 
ternal circuit is open while the poles correspond to the 
short-circuit case. 



Equation (|3C|) shows that the poles of Ytot coincide 
with the zeros of Z (or the poles of the quantum wire 
admittance). It follows from Eqs. (|27]) and ( |l9| ) that 
there is only one set of eigenfrequencies which are com- 
plex with a negative imaginary part corresponding to de- 
caying fluctuations. Several authors have found two sets 
of the eigenfrequencies in the case of short-ijange interact- 
ing electrons for a three-terminal structural or using an 
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other way for the calculation of the observed current .0 
Both sets of eigenfunctions also describe the decaying 
fluctuations. 

Of more interest, in our opinion, are the zeros of 
Ytot(^)- We show that in this case the conditions can 
be found under which the eigenfrequencies are real and 
charge-wave excitations are very slowly decaying. Ac- 
cording to Eq. (30), lt t = when two equations are 
satisfied simultaneously: 



ReZ(n) = , 

i mZ (n) + ^f^\z(n)\ 2 = o. 

1j 



(31) 
(32) 



The first equation is obviously satisfied by the set of real 
frequencies Q n defined by Eq. (23). The second equa- 
tion can also be satisfied at one of frequencies f2„ if the 
capacitance Co is appropriately fitted, for instance, via 
changing the lead areas. The resonant value of Co is 
estimated as 



Co 



2n 2 hvF 



Physically, Eq. (|T|) means the absence of dissipation 
while Eq. ( |32| ) is a resonant condition for the circuit which 
consists of the quantum wire inductance and the inter- 
electrode capacitance. Under this condition, the resonant 
frequency of the charge waves in the wire coincides with 
that of the L c gCo circuit, which results in a strong in- 
crease of the charge- wave amplitude. 



VI. CONCLUSION 

In this paper we have investigated the linear trans- 
port of interacting electrons in a quantum wire of meso- 
scopic length with massive leads. The key point is the 
full enough account of the actual Coulomb interaction 
inside the wire and of three-dimensional electric field in 
the surrounding media. The Coulomb interaction in a 
mesoscopic quantum wire includes both the direct inter- 
action of electrons with each other and their interaction 
via the image charges induced on the leads. We have 
found an exact analytical solution of the problem. This 
has become possible due to (i) use of the bosonization 
technique which is well suited to consider the low-energy 
excitation of ID interacting electrons, and (ii) solving 
the equation of motion for the bosonic phase field by 
expansion in terms of the eigenfunctions of the electron- 
electron interaction operator, which have been found for 
a model case where the electrodes are plates perpendic- 
ular to the wire. 

We have found that the actual Coulomb interaction af- 
fects strongly the electron density distribution along the 
wire in comparison with that in the case of short-range 
interacting electrons in the conventional LL model. The 
nonlocal Coulomb interaction manifests itself first of all 



in the noticeable increase of the charge density in the 
vicinity of the contacts with leads. Here the electron den- 
sity perturbation decreases exponentially with the dis- 
tance from the contact. This effect is essential when the 
frequency is not too high (fi < a^ 1 ^/ 7 !^). 

Another effect of the Coulomb interaction is the renor- 
malization of the charge-wave velocity. Namely, in con- 
trast with the short-range LL model the long-range 
Coulomb interaction causes the frequency-dependent 
renormalization of the charge-wave velocity. This effect 
manifests itself in the frequency dependence of the real 
part of the impedance. With increasing frequency KeZ 
oscillates between two limiting values which are indepen- 
dent of the interaction. These are zero and h/e 2 . The 
fact that KeZ becomes zero is related to the resonances 
of the charge waves along the wire length. Such reso- 
nances occur also when there is no long-range interac- 
tion. The Coulomb interaction causes the resonances 
to be nonequidistant in frequency, which means that 
the charge- wave velocity is frequency-dependent. At low 
frequency the charge-wave velocity is essentially larger 
than the Fermi velocity. With increasing frequency, the 
charge-wave velocity decreases and tends asymptotically 
to the Fermi velocity. 
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APPENDIX A: THE CURRENT BEING 
MEASURED 

The current j mea s being detected by the measuring de- 
vice in the external circuit is equal to the charge flow in 
the leads. Since the currents in the left and right leads 
are obviously equal each other let us consider the current 
in the left electrode. It is equal to the sum of the charge 
flow through the wire j(x = —L/2) and the charge stored 
at the lead surface per unit time: 



jmeas — 3\x=-L/2 



dQ s 
~dt 



(A.l) 



The charge Q s consists of two components. One is the 
external charge caused by the applied voltage if the wire 
is absent, Qext = CoV a , with Co being the mutual capac- 
ity of the electrodes. The other component is the charge 
Qind induced by the charges located within the wire. The 
latter is determined as 



ds a , 
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where 

e dip(r) 
4tt dn ' 

ip is the potential of the charge distributed along the wire 
and n is the outward normal to the lead surface. The po- 
tential if is expressed in terms of the charge density p(r) 
via the Green function G(r,r') of the Laplace equation 
with zero boundary conditions at the lead surfaces. 
Direct calculation results in 



Qind — 



L/2 , / / / \ dGJx,x') 
dx'p{x',t) 

-L/2 



dx 



x = -L/2,q=0 



Using Eq. (Q) and the continuity equation 

dp _ _dj_ 
dt dx ' 

we obtain the induced current 



rfQind ., L 1 



L/2 



L/2 



dx'j(x',t), (A.2) 



where j{x,t) is the particle current in the wire. 
Combining Eqs. (A..1) with (A.2) one obtains 



L 



L/2 dV 

-L/2 ''• 



(A.3) 



The first term in equation ( A.3 ) is a current induced by 
electrons moving in the wire, while th e se cond one is a 



trivial capacitance current. Equation ( A-3 ) is a particu- 
lar case of the general Shockley theorem. For arbitrary 
form and configuration of the leads the Shockley theorem 
is presented as follows 



L /2 

± [ dx'j(x',t)F(x') + C ^, (A A) 

Va J -L/2 at 



where V a is the potential difference between the leads 
and F(x) is the electric field along the electron trajefc 
tory due to V a . According to the original derivation,E3 
the field F(x) appears here as a result of using the reci- 
procity theorem when calculating the charge induced in 
the leads by the charges moving along the wire. Thus, 
F(x) has a sense of the external electric field which does 
not include—the polarization of the ID electron system. 

RecentlyEStJ a question was discussed regarding which 
electric field determines the measured electric current in 
the quantum wire - the external field or the internal one 
- which depends on the polarization of ID electrons. In 
this connection we note t hat i n our case the current cal- 
culated according to Eq. ( |A.4[ ) does not depend on what 
electric field is used. The internal electric field is defined 
by the right-hand side of Eq. (|^) , where the first term is 
the external field F ext and the second one is the induced 



field Fi n d- It is easy to see using Eqs. (p6h, (|l4|), and ( |l5| ) 
that the integral of the product j(x)Fi n d(x) is zero: 



-L/2 



dx 



dx' 



Since in our case F ext is indepe nden t of x, the measured 
current can be found from Eq. ( |A.3| ). 
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